Rank-dependent deactivation in network evolution 
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A rank-dependent deactivation mechanism is introduced to network evolution. The growth dynamics of the 
network is based on a finite memory of individuals, which is implemented by deactivating one site at each time 
step. The model shows striking features of a wide range of real- world networks: power-law degree distribution, 
high clustering coefficient, and disassortative degree correlation. 
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Since 1998, with the small-world model introduced by 
Watts and Strogatz we have witnessed the emergence of 
a new science of networks, which has powerful function of 
describing structures 10, 01 and dynamics [4, 5] of many real 
systems. A network is a mathematical object which consists 
of vertices connected by edges. Despite differences in their 
nature, most real-world networks are characterized by similar 
topological properties, in contrast to those obtained by tra- 
ditional random graphs. For instance, real networks display 
higher clustering than that expected from random networks 
1^. Also, it has been found that many large networks are 
scale free (SF) 120, which means a power-law distribution of 
connectivity, P{k) ^ k^'', where P{k) is the probability that 
a vertex in the network is of degree k and 7 is a positive real 
number determined by the given network. In order to under- 
stand how SF networks arise, Barabasi and Albert (BA) pro- 
posed an evolving network model in 1999, which grows at a 
constant rate and new vertices attach to old ones with proba- 
bility n(fc) ~ fc ITJ. In this way, vertices of high degree are 
more likely to receive further edges from newcomers. 

For most networks, however, aging of sites usually occurs. 
For instance, in reference networks old papers are rarely cited; 
in social networks people of the same age are more likely to 
be friends. To study the effect of aging on network evolution, 
the BA model has been modified by incorporating time depen- 
dence in the network S |1, [M El [il [Tj, Si] . Dorogovtsev 
and Mendes studied the case when the connection probability 
of the new site with an old one is not only proportional to the 
degree k but also to the power of its present age t^" [8J. They 
showed both numerically and analytically that the scale-free 
nature disappears when a < —1. As an alternative, Zhu et al 
introduced the exponential decay function e^^'^ of its present 
age to the BA attachment probability |9i. It was found that 
the produced network is significantly transformed besides the 
change in the degree distribution. On the other hand, Klemm 
and Egufluz 11511 observed the negative coiTelation between 
the vertex age and its rate of acquiring links from the net- 
work of scientific citations. Based on that, they investigated 
the finite collective memory of popular individuals and pro- 
posed a highly clustered scale-free network model lisi |3. 
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The model accounts for three empirical features: preferential 
attachment, power-law degree distribution, and negative cor- 
relation between age and connection rate. 

Recently, Fortunato et al proposed a criterion of network 
growth that explicitly relies on the ranking of vertices ITtIi . 
which originates from the idea that the absolute importance 
(popularity or fitness) of an object is often difficult or impos- 
sible for strangers to measure in social networks. Instead, it 
is quite common to have a clear knowledge about the rela- 
tive values of two objects, i.e., who is more popular or richer 
between two individuals. The rank-driven mechanism gen- 
erates networks with the scale-free degree distribution when 
the probability to link a target vertex is any power-law func- 
tion of its rank, even when one has only partial information of 
vertex ranks iflTlfTsll . Since the perception of how items are 
ranked requires far less information than their actual impor- 
tance, the rank-driven mechanism can well mimic the reality 
in many cases that the relative values of agents are easier to 
access than their absolute values. In this paper, we integrate 
rank with deactivation and study their influences on network 
evolution. Simulations show that interesting statistical proper- 
ties of the generated network display good features observed 
in realistic systems. 

The present model is based on the rank-dependent deacti- 
vating of vertices, which describes the growth dynamics of a 
network with dkected links, run as follows. First, start from 
an initial network of m completely connected seeds, whose 
states are active. By fc'" we denote the in-degree of the vertex, 
i.e., the number of edges pointing to the vertex. At each time 
step, add a new vertex n with m outgoing edges. The new 
vertex is disconnected at first, so fcj" = at this point. Each 
vertex i of the m active vertices receives exactly one incoming 
edge, thereby k'" fcj" + 1. Then activate the new vertex n 
and deactivate one (denoted by j) of the m + 1 active vertices 
with probability 



where a > is a constant bias and 7 is the normalization fac- 
tor Rj £ [1, TO+1] is the rank of j among the rn+l active ver- 
tices. The average connectivity of the network is given by the 
number rn of outgoing edges per vertex. The new added ver- 
tex is always in the active state first and receives edges from 
subsequently added vertices until it is deactivated. Note that 
the larger rank a vertex possesses, the more difficult for it to 
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FIG. 1 : (Color online) Illustration of average R versus average k'" 
of the generated networks for various values of m: m = 10 (solid 
line), 20 (dashed line), 50 (dotted line), and 100 (dashed dotted line), 
respectively. The size of the network is A'' = 10^. 



FIG. 2: (Color online) Log-log plots of the vertex degree distribution 
of the generated networks for different values of a: a = (solid 
line), m (dashed line), and A'^ (dotted line), respectively. All the 
experiment networks have a total number of vertices A'^ — 10^. 



be deactivated. For the case of the citation network, Eq. ([T]l 
means that the famous paper cited mostly is less possibility 
to be forgotten. In Ref. llTll . the model grows according to 
the rank-based preferential attachment n(n ^ i) ~ In 
case vertices are sorted by age R{t) — t, the older the vertex 
is, the higher possibihty for it gaining new edges, coinciding 
with ours. 

The choice of prestige measure can be arbitrary, either topo- 
logical measures or physical ones. In present work we sort 
vertices by age for simplicity. Namely, the older the vertex is, 
the larger rank it possesses. Supposing the distribution of the 
in-degree of active vertices at time t denoted by p(fc'", t), then 
we can write out the master equation 



p{k'" + i,t + !) = [!- n(fc'")]p(fc'", t), 



(2) 



lation between R and A;'". In Fig. [T]we plot i? as a function of 
/c'" of the model by numerical simulations. In order to reduce 
statistical error, the in-degrees of the vertices are calculated as 
an average. As it can be seen, there is a rough power law be- 
tween R and fc'", R ^ ifTgll . Then one easily obtains 
n(A:") - fii^{k'"y-^Il{R), where n(i?) is the deactivation 
probability of a vertex with rank R. Substituting them into 
Eq. we get 



i,t + i) 



fij^{k'Y'^^{R)]p{k'",t) 



p(fc'",<). (3) 



where n(fc™) is the deactivation probability of a vertex with 
in-degree fc'". To do further calculation, we should get the re- 



The subsequent thing is just to follow the analytical method 
in Ref. lilSJ . Imposing the stationarity condition p{k™. t + 



3 



1) =p(k^,t) yields 



+ 1) - Pin = -iM^Tlpi^k'^), 
a + fi[k'^y 

Assuming fc'" changes continuously, we have 



and accordingly obtain the solution 

p(fc'") =6[a + /i(fc'")'^ 



(4) 



(5) 



(6) 



with the appropriate normalization constant h. In case the total 
number N of vertices in the network is large compared with 
the number m of active vertices, the overall in-degree distri- 
bution A^(fc'") can be calculated as the rate of the change of 
the in-degree distribution p(fc") of the active vertices, which 
obeys 



iV(fc" 



dp{k" 



b^fiiy^k" 



[a+^l{k'") 



in\Z''] — 7 — 1 



(7) 



If we chose the value of the bias a — 0, Eq. d?) is no other than 
the probability distribution of the total degree fc = A;'" + m of 
vertices 



N{k) = h-i^r^vk- 



(8) 



In Fig. 12] we plot the total degree distribution of the result- 
ing networks with different values of a. All the plots display 
the good right-skewed behavior, which is reasonably in agree- 
ment with the condition of many realistic systems ll20l l2lll . 
Especially for a = 0, we notice beautiful power laws. 

The "complexity"of networks usually cannot be fully char- 
acterized by the degree distribution of vertices. Instead, the 
self-organization of structures of complex networks is mathe- 
matically encoded in various correlations existing among dif- 
ferent vertices. To describe the network structure in more de- 
tail, several other topological quantities have been introduced 
to the statistics of networks, such as clustering coefficient, de- 
gree correlation, shortest path length, and so on. In the fol- 
lowing, we shall go beyond the degree distribution and discuss 
those quantities. 

Let us start with the clustering coefficient C of the network, 
which is defined as the average probability with which two 
neighbors of a vertex are also neighbors to each other. For 
example, if a vertex i has ki edges, and among its ki nearest 
neighbors there are edges, then the clustering coefficient of 
i is defined by 



2e, 



kiiki 1) 



(9) 



In order to compute the clustering coefficient, we shall con- 
sider the network as undirected and denote hy k — fc'" + m 
the total degree of vertex i. In the deactivation model, new 
edges are created between the active vertices and the added 
one. Moreover, all the active vertices are connected. At each 
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FIG. 3: (Color online) Average clustering coefficient C as a function 
of k (a) and (b), respectively. The parameter m = 50 for all plots. 



time step, the degree ki of each active vertex i increases by 1 
and Ci increase by m — 1. Therefore, the evolutionary dynam- 
ics of ki and ai^e given by 



ki = m + t, 
dci 

—— ~ r 
dt 



1. 



(10) 
(11) 



Integrating Eq. (fTTT i with the boundary condition 6^(0) = 
m{m — l)/2 and substituting the solution into Eq. (|9|l, we 
recover the clustering coefficient C{k) restricted to the ver- 
tices of degree k |[T6ll22ll . 



C(fc) 



2(171 — 1) TO(m — 1) 
k- 1 k{k ~ 1) ' 



(12) 



The expression indicates that the local clustering scales as 
C{k) k^^ for large k. In Fig. Oa) we plot the aver- 
age clustering coefficient C{k) as a function of the vertex 
degree k. The best linear fit gives C(fc) ^ k^^ with expo- 
nent ^ = 0.98(6), which coincides with the prediction of Eq. 
( fT2b . This is the signature of a nontrivial architecture in which 
low-degree vertices generically belong to well interconnected 
communities while high-degree ones are linked to many sites 
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FIG. 4: (Color online) Average nearest-neighbor degree k nn as a 
function of the vertex degree k for different values of the bias a. The 
network parameters are N = 10^ and m = 50. 



that may belong to different groups which are sparsely con- 
nected. In Fig. [3jb) we present three typical curves of the 
average clustering coefficient C versus the network size N. It 
is worth noting that the clustering coefficient of the generated 
networks for all cases is higher than that for the correspond- 
ing one-dimensional regular lattices whose value is 3/4 in the 
limit case i23ll . 

Another commonly studied topological quantity is the de- 
gree correlation (or the mixing pattern), which can be charac- 
terized by analyzing the average degree of nearest neighbors, 
defined by [24,1 



\ij kj . 



(13) 



If knn.i does not show any dependence on the degree of i, 
the network is uncorrelated. In case knn.i is dependent on the 
vertex degree i, there are two kinds of correlations. If knn.i 
increases with k, the network is assortative mixing; i.e., ver- 
tices with high connectivity will connect preferably to highly 
connected ones. If knn.i decreases with k, the network is dis- 
assortative mixing; i.e., vertices with high connectivity will 
connect preferably to lowly connected ones |25]. In Fig. |4]we 
show the simulation results of the average degree of nearest 
neighbors fc„„ as a function of k for different values of a. In 
all cases, the degree correlation in the deactivation model is 
disassortative. 

To complete our study of the model, we finally investigate 
the scaling of the shortest path length and the diameter of the 
network. The shortest path length between two vertices is de- 
fined as the minimum number of intermediate vertices that 
must be traversed to go from vertex to vertex. The average 
shortest path length is the shortest path length averaged over 
all the possible pairs of vertices in the network. On the con- 
trary, the diameter is defined as the largest among the shortest 
paths between any two vertices in the network |22]. In Fig. 
|5] we show the scaling behavior of the average shortest path 
length L and the diameter D of the resulting network. Both 




2000 4000 6000 8000 10000 
N 




2000 4000 6000 8000 10000 

N 



FIG. 5: (Color online) Scaling of the average shortest path length 
(a) and the diameter (b) for different values of the bias a: a — 
(squares), m (closed circles), and TV (triangles), respectively. The 
parameter m — 50 for all plots. 



quantities grow linearly with the network size N similar to 
one-dimensional regular lattices. That is to say, the deactiva- 
tion model does not exhibit small-world properties. It indi- 
cates that there are very few effective shortcuts to reduce the 
distance although the network is highly clustered. However, 
the chainlike structure of the model does not necessarily mean 
large values of L and D. Here L and D are still relatively 
small compared with N, especially for the case a — 0. 

In summary, we have suggested a rank-dependent deactiva- 
tion mechanism and studied its influence on network growth. 
The resulting network shows several good features, (i) The 
degree distribution N{k) of vertices is power law, which in- 
dicates the heterogeneous topology, (ii) The clustering coef- 
ficient is larger than that of one-dimensional regular lattices. 
Besides, the local clustering scales asC(fc) ^ fc^^ for large fc. 
(iii) The decay of the average degree of nearest neighbors fc„„ 
with k characterizes the disassortative mixing pattern, (iv) 
The average shortest path length and the diameter grow with 
the network size, which results from deactivating sites dur- 
ing network evolution. Most of above properties have been 
found very common in realistic systems. We hope that the 
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measurement conducted in present work could be applied to 
real networks in the empirical study. The model we have ex- 
plored, however, is possibly the simplest one in the class of 
rank-dependent deactivation growing networks. One can no- 
tice that in our model the out-degree of the vertices remain 
unchanged during the whole evolving period, which is likely 
unreasonable for citation networks. Furthermore, many net- 
works are intrinsically weighted, their edges having different 



strengths (l^ . Naturally one can generalize the present model 
to the weighted case and sort ranks by the vertex strength or 
fitness. There exists a series of improvements to be made in 
the future. 
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